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ABSTRACT 


A  simple  formula  for  the  Newtonian  drag  coefficient  is  derived 
for  the  blunted  cone-sjdiere  ("ice-cream  cone")  at  small  angles  of  attack.  Graj^s 
are  presented  for  Cp  as  a  function  of  cone  semi-vertex  angle  from  6°  to  30°  and 
"bluntness  ratio"  (nose -radius/base -radius)  ranging  in  tenths  from  0  to  0.  5.  For 
given  values  of  the  body  parameters,  the  wei^t  required  for  the  body  to  attain  a 
prescribed  ballistic  coefficient  is  then  determined.  Fineness  ratios  for  different 
specified  cone -angles  and  degrees  of  bluntness  are  provided  for  reference. 
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Introduction.  Newtonian  impact  theory  makes  two  basic  assumptions: 
first,  that  the  shock  wave  lies  along  tlie  body,  surface;  second,  that  in  striking  die 
body,  the  airstream  loses  the  component  of  momentum  normal  to  the  surface  -  -  the 
flow  is  regarded  as  continuing  along  the  surface  with  the  tangential  component  of 
momentum  unchanged. 

The  local  pressure  coefficient  obtained  under  the  Newtonian  hypothesis 
therefore  depends  only  on  the  local  inclmation  of  the  surface  to  the  airstream: 


C 


P 
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vdiere  and  are  the  free -stream  velocity  and  its  component  normal  to  the  body 

% 

surface  respectively.  Cp  is  taken  to  be  zero  over  those  portions  of  the  surface 
shielded  from  the  flow. 

The  second  of  the  basic  assumptions  made  above  neglects  centrifugal 
forces  associated  with  curvature  of  the  body  surface,  and  the  condition  of  the  first 
assumption  is  attained  only  in  the  limit  as  0.  (M«,  is  the  free-stream 

Mach  number  and  y  the  ratio  of  specific  heats. )  In  this  connection,  the  relevam 
parameter  appears  not  to  be  the  Mach  number  itself  so  much  as  the  ratio  M«,/P, 
where  F  is  the  body  fineness  ratio 

F  -  Total  body>length/Maximum  body- diameter. 


*  See  A]^ndix  I. 
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For  values  of  M^/F  less  than  2  for  example,  Newtonian  estimates  of  the  local 
pressure  coefficient  for  cones  have  been  found  ^  to  be  in  error  by  more  than  17%.  On 
the  other  hand,  Newtonian  theory  has  been  found  to  give  qualitative  and  even  quanti¬ 
tative  information  of  acceptable  accuracy  in  practical  cases  of  interest  at  hypersonic 
speeds,  and  the  following  computations  of  drag  coefficients  and  related  quantities 
should  prove  useful  even  at  higjh  altitudes  when  the  contributions  of  viscous  effects  * 
and  ablation  are  incorporated. 

The  following  sections  are  concerned  with  a  single  configuration  at 
small  angles  of  attack.  The  body  studied  consists  of  a  cone -frustum  capped  at  both 
ends  by  tangent  sjiierical  surfaces.  Cone  semi-vertex  angles  between  6°  and  30°  are 
considered  with  bluntness  ratios  (nose-radius/base-radius)  ranging  in  tenths  from 
0  to  0.  3.  Under  the  assumption  of  negligible  pitching  and  zero  angle  of  attack,  a 
simple  formula  is  derived  for  the  corresponding  Newtonian  drag  coefficients.  This 
was  programmed  for  the  IBM  7090,  and  the  tabulated  values  of  Cq  so  generated  have 
been  presented  in  graphical  form.  Also  presented  are  estimates  (based  on  Cp^)  of  the 
weight  required  for  the  body  to  have  a  prescribed  ballistic  coefficient.  Fineness  ratios 
are  given  for  reference. 

The  Drag  Formula .  The  blunted  cone -sphere  configuration  to  be 
considered  has  rotational  symmetry.  A  coordinate  system  will  be  chosen  with  origin 
at  the  nose  and  x  measured  along  the  axis  of  symmetry.  (See  Figure  1. )  Whh  the 
surface  then  described  by  the  radial  distance  R  =  R  (x),  the  axial -force  coefficient 
may  be  found  by  integration  over  the  surface  actually  exposed  to  the  flow: 


*  These  are  relatively  important  for  sharp-nosed  bodies  during  early  re-entry. 
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Cx  = 


Cp  sin  6  dS 


where  A  is  the  reference -area  and  Q  =  arctan  R'(x).  * 

If  the  angle  of  attack  is  small,  so  that  the  entire  surface  between  the 
planes  x  =  0  and  x  =  -t  is  that  exposed  to  the  flow,  and  if  furthermore  the  angular 
pitch- rate  q  is  negligible  in  comparison  with  the  velocity  v^,  then  integration  and 
simplification  **  lead  to  an  approximate  expression  for  the  drag  coefficient: 


CD-Cx 


R(x)  R’(x)^  ^ 
1  +  R’(x)2 


(0) 


where  the  interval  (o,  is  to  be  broken  into  three  segments  corresponding  to  the 
different  definitions  of  R  over  the  nose,  midsection,  and  base  of  the  body. 


A  spherical  tip  of  radius  R^^  is  generated  by  revolving  a  segment  of  the 
circle  (x  -  R^l  ^  +  [R(x)]^  =  R^j^:  the  radial  distance  for  the  nose  portion  of  the 
body  is  therefore 

R  =  [  2Rn  X  -  x2  ]  for  0  <  X  <  X],  (1) 


*  See  Appendix  1. 

**  See  Appendix  II. 
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where  tangency  occurs  at  Xj.  We  denote  by  17  the  (constant)  semi-vertex  angle  of  the 
conical  midsection,  and  note  that 


xi  =  -  R  sin  T) 

*  n  n  ' 


The  radial  distance  for  the  midsection  is  them  simply 

R  =  [x+xi/sinrj)  tan?)  forx2<x<X2.  (2) 

where  the  spherical  back,  having  radius  Rj^,  is  tangent  at 

X2  =  -t  -  Rjj  sin  7). 

Finally,  the  portion  of  the  spherical  back  which  is  not  shielded  from  the  airstream  has 
radial  distance 

R  =  [Rjj^  -  (x  -  t  for  X2  X  <  -t  ,  (3) 

where  the  maximum  radius  of  the  body  is  Rjj  =  R(-t  ) 

Computation  of  Cp  as  represented  by  equation  (0)  thus  requires  the 
integration  of  three  separate  expressions, R  being  defined  in  the  respective  intervals 
as  noted  above: 
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for  0  s  X  s  xj, 
for  Xj  <  X  <  X21  (4) 

for  ^  X  ^  I  . 

The  second  expression  is  easily  handled  with  R  as  independent  variable  going  between 
Rn  cos  r\  and  Rj,  cos  r\,  and  the  other  two  integrations  are  routine. 

We  shall  find  it  convenient  now  to  take  as  reference-area  the  maximum 
cross-section  of  the  body. 


RR”^ 

1  +R'2 


(Rn  -  x)3/Rj2 


RR'  sin^  7j 


{I  -x)3/Rj^2 


A  =  tRj^. 

Substitution  and  integration  of  (4)  in  equation  (0)  then  yields 

Cjj  =  <2  (1  .  sin4  ^  2(1  -  ic2)  sin^  )?  cos^  Tj  +  sin^  T),  (5) 

where  we  have  introduced  the  "bluntness  ratio"  k  -  Rjj/Rb  as  a  convenient  parameter. 
A  trigonometric  identity  now  gives  a  simple  formula  for  Cq  as  a  function  of  k  and  77 : 

Cq  =  (<2  cos^  77+  [2  -  sin2  77)  sin^  77.  (6) 

The  same  technique  (temonstrated  above  has  been  used  to  find  Cq  for 
the  case  of  a  blunted  cone  which  has  a  flat  back  at  distance  x  »  ^  from  the  iK)se.  The 
corresponding  formula  (obtained  on  replacing  X2  by  in  equation  (4))  is  offered  for 
comparison; 

Cp  =  cos^  q  +  2  sin^  q  (flat  back). 
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The  obvious  similarity  of  the  two  expression  conforms  to  the  expected  similarity  of 
body- shape  when  tj  is  small.  It  should  be  ix>ted  that  the  formula  we  have  derived  for 
the  cone  with  sj^erical  back  differs  from  the  flat-back  case  not  only  in  the  ptresence 
of  a  (sin^  tj)-term  whose  contribution  may  be  negligibly  small  if  K/ri  is  large;  but 
also,  for  a  given  nose -radius,  base-radius,  and  exposed  body-length  t ,  the  spherical 
back  gives  a  larger  semi -cone  angle  i)  than  does  the  flat  back. 

A  final  simplification  of  equation  (6)  gives  a  concise  formula  for  the 
Newtonian  drag  coefficient  of  a  blunted  cone-sphere  at  small  angles  of  attack: 


Cj)  =  I  -  [I  -  ]  COS^  7) 

2 

where  k  =  Rn' »  R-b 


Fineness  Ratio  It  is  easy  to  verify  that  the  semi-vertex  angle  tj  in  a 
blunted  cone -sphere  satisfies  the  relation 

sin  77  =  (Rj^  -  Rjj)/(t  -  Rjj). 

If  the  total  body  length  is  denoted  by  L,  then  t  =  L  -  Rj^,  and  we  may  use  the  fact 
that  the  fineness  ratio  is  F  =  L/2Rj^  to  write  instead 

sin  77  =  (1  -  x)/(2F  -  [fc  +  11), 

determining  the  angle  77  when  the  ratios  of  Rj^,  Rjj,  and  L  are  known.  TTie  relation 
b^een  77  and  F  expressed  by  this  equation  is  presented  graphically  for  different 
values  of  the  bluntness  ratio  k  in  Figure  4. 


It  may  be  of  interest  for  comparison  to  note  that  the  semi -vertex  angle 
T}'  for  a  flat -backed  cone  satisfies  the  relation 

sin  n’  =  (Rb  *=0®  ‘  ^n^’ 

where  for  this  case  t  =  L  It  follows  that  if  a  cone  with  sjdierical  back  were  to  have 
the  same  bluntness  and  fineness  ratios,  it  would  have  to  have  a  (larger)  semi-vertex 
angle  tj,  determined  by  the  ratio 

sin  17 _ 1  -  K _ 

sin  r}‘  ~  (cos  tj'  -  sin  )  -  k 


Weights.  While  nearly  constant  over  considerable  altitude  ranges  in 
the  case  of  blunt-nosed  bodies,  the  ballistic  coefficient  (or  " weight- to»drag"  ratio) 
W/CjyV  may  vary  appreciably  at  hi^  altitudes  for  sharp-nosed  bodies.  Viscous  drag 
is  relatively  insignificant  at  50.000  feet  however,  and  common  usage  seems  to  base 
W/CjjA  on  the  (pressure)  drag  at  this  altitude  when  no  further  qualification  is  mention 
ed.  Since  maximum  deceleration  of  a  high  performance  re-entry  body  might  be  ex¬ 
pected  to  occur  at  an  even  lower  altitude  (determined  by  the  ballistic  coefficient  and 
the  re-entry  path  angle),  the  Newtonian  value  of  Cj-j  may  therefore  prove  useful  in 
estimating  the  weight  required  for  a  specified  ballistic  coefficient  to  be  attained. 

If  the  desired  ballistic  coefficient  is  ^ ,  then  substituting  for  the 
reference -area  A  gives 


W/Rb^  =  IT 
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The  accompanying  grajrfis  of  W/Rjj^  vs  rj  for  different  values  of  k  are  based  for 
scaling  convenience  on  a  ballistic  coefficient  of  100  Ibs/ft^.  Hence,  for  a  given  semi¬ 
vertex  angle  tj  (degrees)^  bluntness  ratio  «,  and  base-radius  R5(ft),  the  weight  W  (lbs) 
necessary  for  a  ballistic  coefficient  /3  (lbs/ft2)  is  given  by 

W  =  . 

where  the  expression  in  brackets  is  the  value  read  off  from  the  appropriate  curve. 


APPENDIX  I 


The  component  of  free -stream  velocity  which  is  in  the  direction 

of  a  normal  h  to  the  surface -element  dS  has  magnitude 


n 


=  cos 
00 


so  the  mass  striking  dS  normally  in  unit  time  is  the  product  of  the  volume  generated 

and  the  free -stream  density  p  : 

00 


dm  =  (Vjj  dS). 

The  rate  of  change  of  momentum  of  the  mass  striking  dS  in  the  direction  of  the 
normal  gives  the  local  pressure -force 

dF  =  (p  -  p^)  dS  =  Vj^dm  =  p„  v^^^  dS  , 

where  p  and  p^  are  the  local  pressure  and  the  free -stream  pressure  respectively. 

As  a  function  of  the  local  pressure  difference  (p  -p^o)  the  free- 
stream  dynamic  pressure 


^  = 


"ttP  v_ 

2  «  * 


the  local  pressure  coefficient 


✓  V 
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thus  depends  only  on  the  local  inclination  of  the  surface  to  the  airstream.  (The 
flow  is  considered  to  be  completely  separated  with  Cp  =  0  over  those  portions  of  the 
surface  shielded  from  the  flow. ) 

We  are  concerned  with  a  smooth  body  having  axial  symmetry.  Its 
surface,  accordingly,  may  be  described  by  giving  the  radial  distance  R  as  a  differ¬ 
entiable  function  of  just  the  distance  x  measured  from  the  nose  of  the  body  (as  origin) 
along  the  axis  of  symmetry  The  local  inclination  of  the  surface  with  respect  to 
this  axis  is  then 


0  =  0(x)  =  arctanR'(x) 

When  there  is  an  angle  of  attack  a,  the  local  pressure  coefficient 
varies  not  only  with  x  but  also  with  the  angle  w  measured  in  a  plane  orthogonal  to 
the  axis.  (See  Figure  I )  Cp  is  thus  dependent  on  a,  d,  and  u.  * 

The  above  relations  being  understood,  the  force  acting  on  an  area- 
element  dS  may  be  expressed  simply  as 

dF  =  q,  Cp  dS 

The  component  of  this  force  in  the  direction  of  the  axis  is  then  (sin  0)dF,  and 
integrating  over  the  portion  of  the  surface  actually  exposed  to  the  airstream  gives 
the  total  axial  force 


Cp  sin  0  dS. 


*  See  Appendix  11. 
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The  non-dimensional  axial -force  coefficient  defined  with  respect  to  the 
reference -area  A  by  the  relation 


may  therefore  be  represented  in  the  followmg  form: 


C 


X 


sin  6  dS. 


Expressing  the  element  of  area  dS  by  (  Rdw)  (dx/cos0)  and  substituting  formally 
then  gives 


C 


X 


Cp  RR’dw^dx 


with  integration  limits  corresponding  to  the  surface  exposed  to  the  flow. 
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APPENDIX  II 


The  longitudinal  axis  of  symmetry  of  the  body,  directed  positively 
back  from  the  nose,  has  been  chosen  as  the  X-axis  of  the  body  coordinate  system. 
The  free -stream  velocity  vector  is  assumed  to  be  at  a  small  angle  of  attack  a  (of 
no  greater  magnitude  than  ij)  relative  to  the  X-axis,  determining  with  it  a  plane  in 
which  pitching  occurs  about  the  center  of  gravity  x^g  at  an  angular  rate  q. 

With  these  conventions,  the  local  normal  component  of  velocity 
Vn  may  be  expressed^’  ^  in  the  form 


V  =  V  cos  a  [  sin  0  -  tan  a  sin  u  cos  0  ]  -  q  sin  w  [  (x  -  x„g)  cos  0  +  R  sin  0  ] , 
n  o 


or  equivalently,  upon  substituting  for  cos  0  and  tan  0, 


,  ,  V«,  cos  O' 

V  (x,  w)  =  ,  , '  , 

”  ^  i+r-2  L 


R'  -  sin 


- * .  (  R R'  +  X  -  x-J  +  tan 

Voo  cos  a  ‘  eg' 


It  then  follows  from  definition  of  Cp  that 


Cp  (X,  a,  w)  =  --  (R'  -  G  sin 

1  +  R’2 


where 


_ 

v^  cos  or 


G(x) 


1  R  R’  +  X  -  x^g]  +  tan  a. 
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The  entire  surface  corresponding  to  0  <  x  <  -t  and  0  s  oi  ^  2jr  is 
exposed  to  the  flow  at  small  angles  of  attack  (a  <  rj),  so  the  surface  integral  defining 
Cjj  may  be  evaluated  as  an  iterated  integral  with  these  limits  on  x  and  w.  Taking 
advantage  of  symmetry,  we  obtain 

2  r 

Cy  =  T“  \  \  R  (x)  R'  (x)  C  (x,  a,  u)  du  dx. 

^  ^  0  -ff/2  P 


It 


remains  to  substitute  for  C 


P 


in  the  form  obtained  above: 


>+ir/2 
-  jr/2 


(R‘  -Gsin  u;)^  du)  =  (2R*^  +  G^) 


'  4ir  C^'  RR’^  .1  2  .  2ir  r^RR'  ,2 

"  Cx  =  T  \  -  ^  cos^  “  T"  \  - 7 

^  ^  *^0  1  +  r'2  a  ^  0  1  +R'2 


dx. 


The  drag-force  is  directly  opposed  to  the  body’s  velocity -vector  (relative  to  the  flow) 

at  an  angle  a  from  the  axis,  and  the  integrand  on  the  rij^t  is  clearly  bounded  with  q. 

2 

Hence  the  relation  cos  a  *»  1  for  small  a  and  the  assumption  that  (q/v^)  is  negli¬ 
gible  leaves  just  the  bracketted  expression  asserted  in  equation  (0) 


3-22-3790 
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a 


X,  R,  0) 

R(x) 

R'  (X) 


*cg 

R„>  Rv 


I 

L 

F 


V 

n 


A 


List  of  Symbols 

angle  of  attack  (front  body  axis  to  relative -velocity  vector) 
cylindrical  coordinates,  x  along  axis  of  symmetry 
body -radius  in  (R,  w)  -  plane 

rate  of  change  of  body-radius  along  axis  of  symmetry 
body  center  of  gravity  (distance  from  nose) 
nose-radius,  base-radius  of  body 

"bluntness  ratio",  x  =  R  /R^ 

n  D 

total  length  exposed  to  flow,  "wetted"  length 
total  length  of  body 
fineness  ratio,  F  =  L/2Rj^ 

semi -vertex  angle  for  conical  midsection 
local  inclination  of  body  surface,  d  =  arctan  R’  (x) 
angular  pitch- rate  about  transverse  axis  through  Xj,g 
local  pressure,  local  pressure  coefficient 

free- stream  pressure,  free -stream  density 

free -stream  velocity,  free  stream  Mach  number 

dynamic  pressure,  ^  Poo'^,» 

component  of  free-stream  velocity  normal  to  body  surface 

reference -area  for  non-dimensional  coefficients 
axial-force  coefficient,  Cj^  =  axial  force/q^A 

drag-force  coefficient,  Cj^  =  drag  force/q^  A 
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Distribution  List  (continued) 

Group  42 

W.  W  Ward 
P  A  Northrop 
C  B  Slade 
F  Betts 
J  B  Cooper 
j  B  Jelatis 
B  H  Labitt 
V  C  Martins 
W  R  Renhult 
R  G  Sandhcilm 
A  F  Standing 
C  M  Steinmetz 
R  C  Yost 

Group  43 

J  G  Biirry 

Group  46 

C  W  Jones 

Group  47 


D  L  Clark 

Group  4S 

R  C  Butman 

Group  7 1 

F  G.  DeSantis 
P  Alberti 

Group  744 


J  F  Hutzenlaub 


